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Abstract 

We consider the problem of enumerating spanning trees on lattices. Closed- form 
expressions are obtained for the spanning tree generating function for a hypercubic 
lattice of size A'^i x x • • • x N4 in d dimensions under free, periodic, and a combination 
of free and periodic boundary conditions. Results are also obtained for a simple quartic 
net embedded on two non-orientable surfaces, a Mobius strip and the Klein bottle. 
Our results are based on the use of a formula expressing the spanning tree generating 
function in terms of the eigenvalues of an associated tree matrix. An elementary 
derivation of this formula is given. 
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1 INTRODUCTION 

The problem of enumerating spanning trees on a graph was first considered by Kirchhoff [0 
in his analysis of electrical networks. Consider a graph G = {V, E} consisting of a vertex set 
V and an edge set E. We shall assume that G is connected. A subset of edges T C -E is a 
spanning tree if it has — 1 edges with at least one edge incident at each vertex. Therefore 
T has no cycles. In ensuing discussions we shall use T to also denote the spanning tree. 

Number the vertices from 1 to |V^| and associate to the edge Cij connecting vertices i and 
j a weight Xij, with the convention of xu = 0. The enumeration of spanning trees concerns 
with the evaluation of the tree generating function 

T{G; {x,,}) = E n (1) 

TCE CjjGT 

where the summation is taken over all spanning trees T. Particularly, the number of spanning 
trees on G is obtained by setting 

Nspt{G)=T{G;1). (2) 

Considerations of spanning tree also arise in statistical physics in the enumeration of 
close-packed dimers (perfect matchings) . Using a similar consideration, for example, one 



1 



of us 0] has evaluated the number of spanning trees for the simple quartic, triangular and 
honeycomb lattices in the limit of |\^| — > oo. In this Letter we report new results on the 
evaluation of the generating function Eq. (|l|) for finite hypercubic lattices in arbitrary di- 
mensions. Results are also obtained for a simple quartic net embedded on two non-orientable 
surfaces, the Mobius strip and the Klein bottle. As the main formula used in this Letter is a 
relation expressing the tree generating function in terms of the eigenvalues of an associated 
tree matrix, for completeness we give an elementary derivation of this formula. 

2 THE TREE MATRIX 

For a given graph G = {V, E} consider a \V\ x \V\ matrix M(G) with elements 

[ E^k=i^ik, « = J = 1,2,---,|V"| 
Mij{G) = < —Xij, if vertices i, j, i ^ j, are connected by an edge (3) 

[ 0, otherwise. 

We shall refer to M(G') simply as the tree matrix. It is well-known 0, ^ that the tree 
generating function, Eq. (|lD, is given by the cofactor of any element of the tree matrix, and 
that the cofactor is the same for all elements. Namely, we have the identity 

T{G; {xij}) = the cofactor of any element of the matrix M(G). (4) 

The tree generating function can also be expressed in terms of the eigenvalues of the 
tree matrix M(G'). We give here an elementary derivation of this result which we use in 
subsequent sections. 

Let M(G') be the tree matrix of a graph G = {V, E}. Since the sum of all elements in a 
row of M(G) equals to zero, M(G) has as an eigenvalue and, by definition, we have 

det \Mij{G) - X6,j\ = -XF{X) (5) 

where 

F{X) = lli^.-X), (6) 

1=2 

A2, A3, ... , A|y| being the remaining eigenvalues. 

Now the sum of all elements in a row of the determinant \Mij{G) — XSij\ is —A. This per- 
mits us to replace the first column of det\Mij{G) — X6ij \ by a column of elements —A without 
affecting its value. Next we carry out a Laplace expansion of the resulting determinant along 
the modified column, obtaining 

\y\ 

det\M,^{G) - X6,,\ = -A^Ca(A), (7) 

1=1 

where Cji(A) is the cofactor of the {il)-th element of the determinant. Combining Eqs. 
(S)"(0)) led to the identity 

FiX)=Y:Ca{X). (8) 

i=l 



2 



Now, Cji(O) is precisely the cofactor of the (il)-th element of M(G') which, by Eq. (1), is 
equal to the tree generating function T(G; {xij}). It follows that, after setting A = in Eq. 
(1), we obtain an expression giving the tree generating function in terms of the eigenvalues 
of the tree matrix [7, p. 39] 

T{G;{x,,}) = —l[K. (9) 

This result can also be deduced by considering the tree matrix of a graph obtained from G 
by adding an auxiliary vertex connected to all vertices with edges of weight x, followed by 
taking the limit of x — + [Q. 



3 HYPERCUBIC LATTICES 

We now deduce the closed-form expression for the tree generating function for a hypercubic 
lattice in d dimensions under various boundary conditions. 

3.1. Free boundary conditions 

THEOREM 1. Let be a d- dimensional hypercubic lattice of size Ni x N2 x ■ ■ ■ x with 
edge weights along the ith direction, i = 1,2, . . . ,d. The tree generating function for is 



n n 



AT 



ni=0 



"d=0 



i=l 



^ Xj ( 1 - COS 



(10) 



(ni,...,nrf) (0,...,0), 



where M = NiN2---Nd. 
PROOF. 

The tree matrix of Z^ assumes the form of a linear combination of direct products of 
smaller matrices, 



i=l 

where Jat is an x identity matrix and Hn is the N x N tri-diagonal matrix 



In. 



H 



/ 1 1 
10 10 
10 1 



N 





V 



\ 





1 1 
11/ 



It is readily verified that Hn is diagonalized by the similarity transformation 

SnHnSn^ = An, 



fill 



:i2) 



(13) 
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where Sn and S^^ are N x N matrices with elements 



N) 



cos 



, f'mTx\ 



+ 



m, n = 0, 1 — 1, 



and Ajv is an X diagonal matrix with diagonal elements 



An = 2 cos— , n = 0, 1, . . . , - 1. 



(14) 



(15) 



Here 5m,n is the Kronecker delta. It follows that M(Zd) is diagonalized by the similarity 
transformation 

(16) 
(17) 



S^M(Z,)S^i = Aat, 



where 



and A_v is an AT x A/" diagonal matrix with diagonal elements 



^ni,...,na 

Now, we have Ar 



d 
i=l 



COS ■ 



0,1 



, X , . . . , Z Y J 



A^,; - 1. 



:i^ 



^ni,...,n^ 



for rii = n2 



• • = = 0. This establishes Theorem 1 after 
using Eq. (|). Q.E.D. 

REMARK. The result Eq. ( |7^ generalizes the d = 2 eigenvalues of M(Z2) for Xi = 1 
reported in [7, p. 74]. 

3.2. Periodic boundary conditions 

In applications in physics one often requires periodic boundary conditions depicted by 
the condition that two "boundary" vertices at coordinates (. . . ,?T,j = 1, . . .) and (. . . ,nj = 
A^j, . . .), i = 1,2, . . . ,d, are connected by an extra edge. This leads to a lattice Z^®'' which is 
a regular graph with degree 2d at all vertices. For d = 2, for example, Z^^"^ can be regarded 
as being embedded on the surface of a torus. 

THEOREM 2. Let TF^^ he a hypercubic lattice in d dimensions of size Ni x N2 x ■ ■ ■ x 
with edge weights Xi along the ith direction, i = 1,2, ...,d with periodic boundary conditions. 
The tree generating function for Z^^^ is 



r(Zr ; {x.}) 



^Af^l Ni-l Nd-1 

n n 

ni=0 'T.d=0 Li=l 



AT 



^ Xi ( 1 - cos 



2r2j7r 



rid) 7^ (0,...,0). 



(19) 



PROOF. 

The tree matrix assumes the form 



MfZ 



Per\ 



J2xi[2lN,®l 



N2 



i=l 



/iV,_i ® Gn, ® /7V,+i 



(20) 
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where is the N x N cychc matrix 



G 



/ 1 
10 10 
10 1 



N 





V 1 



1 \ 





1 1 

1 ; 



(21) 



As in Eq. (|T6l), the matrix M(Z^"^'') can be diagonalized by a similarity transformation 
generated by 

Rat = Rni ® Rn2 ® • ■ ■ ® Rnj, 



where Rm is a.n N x N matrix with elements 



where * denotes the complex conjugate, yielding eigenvalues of Gat as 

2mT 



X„ = 2 cos ■ 



N 



n = 0,1 N-1. 



(22) 
(23) 

(24) 
Q.E.D. 



This establishes Theorem 2 after using Eq. (^). 

3. 3. Periodic boundary conditions along m < d directions 

COROLLARY. Let zj^'^"^'-'"'' be a hypercubic lattice in d dimensions of size NiX N2X ■ ■ ■ x 
with periodic boundary conditions in directions 1,2, ... ,m < d and free boundaries in the 
remaining d — m directions. The tree generating function is 



T(Zr^™);{x.}) 



n n 



i=l 



^ ( 1 - COS 



2nj7r 



+ E a;Wl-cos — 

i=m+l ^ 



ni=0 n^=0 

(ni,...,nrf)^(0,...,0). 



(25) 



4 THE MOBIUS STRIP AND THE KLEIN BOTTLE 

Due to the interplay with the conformal field theory H| , it is of current interest in statistical 



physics to study lattice systems on non-orientable surfaces [|T0|, |ri[]. Here, we consider two 
such surfaces, the Mobius strip and the Klein bottle, and obtain the respective tree generating 
functions. 

4.1. The Mobius strip 

THEOREM 3. Let Zf"^ be an M x N simple quartic net embedded on a Mobius strip 
forming a Mobius net of width M and twisted in the direction N , with edge weights Xi and 
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X2 along directions M and N respectively. The tree generating function for Z 



Mob 



IS 



T(Z^°^{xi,X2}) 



n n 



— a;2 1 — cos ■ 



An -3- (-1)'' 



mn 

a:i(l-cos — 



2N 



-ir 



{m,n) ^ (0,0). 



(26) 



PROOF. 

Specifically, let the the two vertices at coordinates {m, 1} and {M— m, A^}, m = 1, 2, ■ ■ ■ , M 
be connected with a lattice edge of weight X2- Then the tree matrix assumes the form 



M(Z 



Mob\ 



2(xi + X'2)Im ®In- xiHm (S) In - X2[Im ® Fn + Jm ® 



N 



(27) 



where 



/ 1 
10 10 
10 1 





V 



\ 





J, 



/ 





M 



K 



1 1 

1 y 

/ 







N 



1 

VI 

1 \ 





1 \ 

1 

1 





; 



■■■ 

V 1 ■■■ / 

Since Hm and Jm commute, they can be simultaneously diagonalized by applying the simi- 
larity transformation Eq. (|l3|) . The transformed matrix SArM(Z^°^)Sj^^ is "block diagonal" 
with N X N blocks 



2 ( xi - xi cos — + X2]In 



X2 {Fj 



N 



m 



0,1 M- 1. 



(2J 



Now, the eigenvalues of Gjy = Fjy + Kjq and F^ — Kjq are, respectively, 2cos[2(n + \)ti/N] 
and 2cos[(2n + l)7r/A^],n = 0, 1,...,A^ — 1. Theorem 3 is established by combining these 
results with Eq. (1). Q.E.D. 

REMARK. For M = 2 and xi = X2 = 1, Eq. gives the number of spanning trees on a 
2 X N Mobius ladder as 



2N-1 



SPT 



2N n 

r 



-r 



2 cos — 

A^ 



2 + (2 + VIY + (2 - V3) 



N 



(29) 



These two equivalent expressions have previously been given by [7, p. 218] and by Guy and 
Harary JT^/ , respectively. 
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4.2. The Klem bottle 

The embedding of an M x simple quartic net on a Klein bottle is accomplished by 
further imposing a periodic boundary condition to Z^°^ in the M direction, namely, by 
connecting vertices of Z^°^ at coordinates {1, n} and {M, n}, n — 1,2, . . . , N with an edge 
of weight xi. This leads to a lattice Zf ^™ of the topology of a Klein bottle. 

THEOREM 4 The tree generating function for Zf-'^™ (described in the above) is 



T{Zf'^'-;{x,,X2}) 



tMN-l 



MN 



n (l - cos ^ j 



X n n 

m=l n=0 



2xi ( 1 — cos 



2m7r \ / riTT \ 

+ zxo 1 — cos — — 



M 



X 



n^-o' [4^1 - 2z, (l 



cos 



(2n+l) 



N 



for M even 
for M odd. 



(30) 



where [n] is the integral part of n. 
PROOF. 



The tree matrix of Z^'^''^ assumes the form 



M(Z^^^''') = 2{xi + X2)Im ®In- XiGm ^ In - X2[Im ®Fj^ + Jm®K_ 



(31) 



To obtain its eigenvalues, we first apply the similarity transformation generated by Rm 
in the M subspace. While this diagonalizes Gm with eigenvalues 2 cos(2m7r/M), m = 
0,1,...,M — 1, it transforms the matrix Jm into 







/ 1 










... 


\ 

















... 


uj 

















... 





RmJmRm 




























-2 


... 









lo 









... 


0/ 


where u = ^'^^/^ ^ and thus M(Z|^i^'°) 


into 












( Ao + Bo 


















\ 
































A2 







B2 















Bm-2 







Am-2 







I 


Bm-1 













Am-\ ) 



(32) 



(33) 



where and are N x N matrices given by 



Am — 2 



xi + X2 — xi cos 



2mTT 



Br, 



-u;"'x2Km, 



M 
m — 0,1, 



N 



X2FN, 

,M- 1. 



(34) 
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The matrix Eq. ( p3D is block diagonal with blocks 



U^{0) = Ao + Bo = X2{2I^-Gr,) 

M -1 



U2N{m) = I J^"" ), m = l,2, 



(35) 



and for M = even, 



Un{M/2) = Am/2 + Bm/2, 

= 2{2xi + X2)In-X2{Fn-Kn), (36) 

where the subscripts of the U matrices denote the matrix dimensions. It follows that we 
need only to find the eigenvalues of the U matrices. 

Eigenvalues of Un{S^) and Un{M/2) can be deduced from those of Gn and Fn — Kn- Fur- 
thermore, eigenvalues of U2N{fn) are obtained from those of On after applying the similarity 
transformation 

/ 2m7r\ 

T2Nirn)U2Nim)T2N{m) = 2 xi + X2 - Xi cos I2N - X2G2N (37) 

V m / 



where 



T2N{m)={ ,,lr\- (38) 



Combining these results with Eq. (||), we are led to Eq. (0) and the theorem. Q.E.D. 
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